Ordered clusters and dynamical states of particles in a vibrated fluid 



O 

o 

(N 



o 



i 

c 

o 
o 



> 

(N 

(N 
O 
(N 
O 



c 

o 
o 



13 



Greg A. Voth 1 , B. Bigger 1 , M.R. Buckley 1 , W. Losert 1 '* , M.P. Brenner 2 , H.A. Stone 2 , and J.P. Gollub 

1 Department of Physics, Haverford College, Haverford PA 1904.1, U.S.A. 
2 Division of Applied Sciences, Harvard University, Cambridge MA 02138, U.S. A 
3 Department of Physics, University of Pennsylvania, Philadelphia PA 19104, U.S. A 
* Present address: IREAP, University of Maryland, College Park, MD 20742-3511 

(February 6, 2008) 

Fluid-mediated interactions between particles in a vibrating fluid lead to both long range attraction 
and short range repulsion. The resulting patterns include hexagonally ordered micro-crystallites, 
time-periodic structures, and chaotic fluctuating patterns with complex dynamics. A model based 
on streaming flow gives a good quantitative account of the attractive part of the interaction. 
PACS: 45.70.Qj,47.15.Cb, 05.65.+b 
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Many important systems consist of granular material 
flowing while immersed in fluid. In early experiments, 
Bagnold used fluid saturated granular flows to remove 
the effects of gravity Q . Others have avoided interstitial 
fluid (even removing air) in order to observe pure granu- 
lar behavior. However, the fluid can produce interactions 
between particles resulting in interesting macroscopic ef- 
fects. For example, granular heaping due to fluid inter- 
actions has been investigated in vertically vibrated deep 
granular beds Also, a variety of interaction effects 

are important in sedimentation (for example |I| ) . Several 
systems show combinations of attractive and repulsive in- 
teractions which result in ordering of particles |^-||] • 

In this letter, we present experimental studies of novel 
phenomena that occur as a consequence of both attrac- 
tive and repulsive interactions between non-brownian 
particles when they are vibrated in a fluid. The obser- 
vations include clustering, ordered crystalline patterns, 
and dynamical fluctuating states. The fluid mediated in- 
teractions can be tuned to produce a striking variety of 
dynamical phenomena, and this provides an interesting 
new mechanism for self-assembly of ordered structures. 

The attractive interaction can be understood quanti- 
tatively using a theory based on the mean streaming flow 
generated by the oscillating particles. A short-range re- 
pulsive interaction is also observed at large vibrational 
acceleration. This repulsion shows a very sharp onset 
as the acceleration is increased. Together the combina- 
tion of attraction and repulsion results in particles being 
bound together without contact over a range of parame- 
ters. We demonstrate how these interactions allow micro- 
crystallites to form, for example hexagons surrounding a 
particle center. Small numbers of particles can form sta- 
ble structures, while larger numbers of particles (> 7) 
move chaotically in a bound state with no long range 
order, a "mesoscopic liquid" . 

The experimental setup in shown in Fig. 1. We con- 
duct all experiments in a rigid 6 cm diameter by 1.5 
cm tall cylindrical aluminum container vibrated verti- 
cally by an electromagnetic vibrator. We completely fill 
the container with a water/glycerol mixture of kinematic 



viscosity v = 8 cS (density =1.1 5 gm/cm 3 ) and a sub- 
monolayer of uniform stainless steel spheres of diameter 
d = 0.794 mm and density 8.0 gm/cm 3 . A glass window 
seals the top of the cylinder to provide optical access and 
to avoid surface waves. The bottom plate is made of glass 
since its smoothness reduces random horizontal particle 
motion. We have used both flat glass and a concave lens 
for the bottom surface. The concave lens (focal length -1 
m) makes no measurable change in the particle interac- 
tions; it is used to keep particles from slowly drifting to 
the edges of the cell due to imperfect leveling or slightly 
non-linear vibrational motion. 

Both the frequency (u> — 2%f) and amplitude (M) of 
the vibration of the container are controlled externally. 
The nondimensional acceleration is given by T = aj 2 M/ g. 
The container oscillations cause the particles to vibrate 
vertically, typically contacting the cell bottom once each 
cycle. The particles remain within 1 or 2 particle di- 
ameters from the bottom over a wide range of acceler- 
ations. We image the entire system from above using a 
fast CCD camera. When the system is illuminated from 
an oblique angle, it is possible to measure the positions 
of the shadows of particles, and thereby to determine the 
amplitude of their vertical motion, A, defined as half the 
peak to peak displacement. The particle Reynolds num- 
ber, Re = dAuj/v, ranges from 2 to 10. 

In Fig. 2 we show the time evolution of an initially 
random distribution of particles. After the vibrator is 
started, the particles quickly collect into localized clus- 
ters, and the clusters then slowly coalesce in a manner 
reminiscent of coarsening in phase transitions. 

A two particle system provides a particularly simple 
flow in which the particle attraction can be compared 
with theory. Figure 3 shows the distance between the 
centers of two particles as a function of time. The ap- 
proach rate increases until they come into contact. The 
shape of these curves provides clear evidence that the 
clustering is due to a fluid-mediated interaction, and not 
to inelastic collisions, which cause clustering in a similar 
system with no fluid 0. The separation also shows a 
small periodic modulation with the period of the vertical 
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vibration. Particles separate near their lowest point and 
approach near the highest point of their motion. 

The attractive interaction can be quantified by consid- 
ering the flow produced by each individual particle's mo- 
tion. The vibrations produce an oscillatory flow around 
each particle, which is well approximated as a potential 
flow, outside of a boundary layer with approximate thick- 
ness S osc = yj2v/u). In the experiment at / = 50 Hz, 
Sosc = 0.2 mm, which is smaller than the particle radius, 
a = 0.397 mm. The potential flow near each oscillating 
particle is u = V</>, where <f> = — \Auj s'm(ujt)a 3 cos(0)/r 2 ; 
here the origin is at the middle of the particle, and 
8 is the angle between a given location on the sphere 
and the forcing direction. At the particle surface, the 
component of the velocity field parallel to its surface is 
it||(0) = \Ausin(wt)sin{6) 0. 

Long ago Lord Rayleigh pointed out that when the 
magnitude of the oscillatory flow u\\ = u\\(6) varies along 
a solid surface, such as that of the particle, a steady sec- 
ondary flow is generated pj]Jl^] . For flow near the sur- 
face, mass conservation then requires that there is also 
a flow perpendicular to the boundary with magnitude 
Mi ~ 5 osc d\\u\\. Since this flow is not entirely out of 
phase with u\\ in the boundary layer, every oscillation 
cycle transports a finite amount of momentum into the 
boundary layer. Hence there is a time independent force 
density on the fluid in the boundary layer, which is par- 
allel to the boundary and of order ^(u||9||U||), where p 
is the fluid density and (•) denotes time average over a 
cycle. This forcing produces a steady flow, which when 
balanced against the viscous force pvu s t ea dy/5 2 sc gives 
the magnitude of the steady flow at the edge of the os- 
cillatory boundary layer u s t ea dy ~ u \\d\\u\\/u>. 

Hence, a steady flow is produced with a magnitude 
Usteady ~ A 2 uo/a. A careful analysis [[l3| shows that this 
steady flow pushes fluid away from the poles of each par- 
ticle; and so there is a perpendicular inflow velocity to- 
wards the equator. This steady inflow is the origin of the 
attractive interactions between two particles. 
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FIG. 1. Schematic of the apparatus in which particles are 
vibrated vertically in a viscous fluid. The vertical height of 
the cell is 19 times the particle diameter, so it approximates 
particles bouncing in a semi-infinite domain. 




El ELI 

' JSP* w» 
* 




FIG. 2. Time evolution of an initially random distribution 
of beads. The attraction created by the streaming flow quickly 
collects the particles into clusters, (f—50 Hz, T=4.5) a) t=0 
s; b) t=8 s, c) t=16 s, d) t=32 s. 



To compute the rate at which particles come together, 
we must determine the inflow velocity far from the parti- 
cle. Besides the oscillatory boundary layer, there is also a 
boundary layer caused by the steady flow itself Q , with 
a scale S stea dy = a/^Resteady, where Re stC ady = A 2 lo/v 
is the Reynolds number of the steady secondary flow. A 
matching argument, connecting the flows in the bound- 
ary layers to a potential flow far away, determines the 
inflow velocity to be v(r) = — 0.53A y /uJi'a 2 /r 3 at a dis- 
tance r from the center of the sphere. Thus, if R(t) de- 
notes the distance between the particle centers, it follows 
that dR/dt = 2v(R) since the particles follow the slow 
changes in the horizontal motion of the fluid. This im- 
plies that the separation between the two spheres should 
decrease according to the law 



R(t) = {R% - 4.2T4 V^a 2 t) 1/4 . 



(1) 



This formula assumes that (i) S osc /S ste ady < F which can 
be rewritten A < a/>/2; (ii) the particles are far enough 
apart that they do not affect each other's boundary lay- 
ers; and (iii) neglects the influence of the bottom plate. 
Our estimates indicate that assumption (i) holds when 
F is small or / is large; (ii) holds except at the final 
moments of approach; and (iii) holds except during the 
small interval in each cycle when the particle bounces off 
the plate. 

In Fig. 3, the central dashed line shows the approach 
curves predicted by Eq. |l| for the conditions of the experi- 
ment. The agreement of the data with both the predicted 
functional form and rate of approach is quite good con- 
sidering that all the parameters used by the theory are 
independently determined. At 20 Hz (Fig. 3a), the mea- 
sured trajectories approach somewhat more slowly than 
the theory predicts. At 50 Hz (Fig 3b), the measured ap- 
proach rate is closer to the prediction, and the functional 
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FIG. 3. Distance between two particles as they are brought 
together by the attractive flow. The solid curves are differ- 
ent experimental runs, and the central dashed curve shows 
the approach predicted by the steady streaming theory. Up- 
per and lower dashed curves show the effect of measurement 
uncertainty in A. The time axes have been shifted so that 
the separation extrapolates to zero at t = 0. (a) f—20 Hz, 
F = 2.9, A = 0.41 ± 0.04 mm. Here the vertical particle mo- 
tion is periodic, (b) /=50 Hz, F = 4.6, A = 0.15 ± 0.04 mm. 
The vertical motion is now chaotic and the mean oscillation 
amplitude is used for A in plotting the theoretical curve. 



form agrees better, most likely because the oscillatory 
boundary layer is thinner in this case and so assumption 
(i) above is better satisfied. The upper and lower dashed 
curves show the effect of measurment uncertainty in A. 
Since at 50 Hz the amplitude is very small and the ver- 
tical motion is chaotic, the uncertainty is rather large in 
this case. Overall, the theory quite accurately captures 
the approach of two particles, and it seems unambiguous 
that the attraction is caused by the streaming mecha- 
nism. 

At larger accelerations, a short-range repulsion be- 
comes prominent and limits the closeness of approach. 
This can be seen in Fig. 4, which displays the steady 
state separation between two particles as a function of 
r for several different forcing frequencies. At each fre- 
quency, there is a step-like rise in separation at a char- 
acteristic r that moves to higher T as the frequency in- 
creases. It is interesting to note (Figure 4b) that there is 
no sharp change in the amplitude of the vertical motion 
of the particles with increasing T, so the onset of non- 
zero separation must be due to the fluid flow and not to 
changes in the bouncing dynamics. The onset of separa- 
tion at each frequency appears when the particles have a 
peak-to-peak vertical amplitude of approximately I Ad. 

A fully quantitative theory of this effect is not avail- 
able. One possible explanation is that the repulsion be- 
comes important when the oscillatory boundary layer is 
thicker than the steady boundary layer. An alternate 
possibility suggested by flow visualization is that the re- 




FIG. 4. (a) For two particles, the equilibrium distance be- 
tween particle centers as a function of acceleration shows a 
rapid transition from contact to a nonzero separation. Sep- 
aration is normalized by the particle diameter, d. (b) Mea- 
surements of the height that particles bounce above the plate 
show no dramatic change at the transition. Frequencies are 
17 Hz (open circles), 20 Hz (filled squares), and 23 Hz (open 
triangles). Larger frequencies have the transition at larger 
acceleration. 



pulsion is due to recirculating zones near the particles cre- 
ated by deflection of the downward part of the streaming 
flow by the plate. Explaining the details of the repulsive 
interaction remains as a challenge for future work. 

Systems of more than two particles display a range 
of fascinating patterns and dynamics, including ordered 
crystallites, time-dependent ordered patterns, and bound 
states with complex particle motion. For example, one 
might expect that when attraction and repulsion are both 
important, the three particle system would form a stable 
symmetric triangular configuration. For low acceleration 
this is nearly the case, but as the acceleration increases, 
the system breaks symmetry to a state with two paired 
and one distant particle as shown in Fig. 5(a). This state 
is fairly stable, but over long times the distant particle 
can wander close enough to cause a new pairing. 

In a system with 7 particles, stable hexagonal states 
form over the range 2.8 < T < 3.0 (Fig. 5b). At slightly 
higher acceleration, the system switches to a time depen- 
dent state with 2 central particles rotating inside a ring 
of 5 others (Fig 5c) . Particle motion here is periodic and 
quite stable, and an animation of this "dance" is avail- 
able online 15 1. At yet higher acceleration, 7 particle 
clusters exhibit apparently chaotic motion. 

For larger numbers of particles, the presence of both at- 
traction and repulsion leads to even more complex many 
body effects, which are generally time-dependent. Fig- 
ures 5(d-f) show a system of many particles at 3 different 
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FIG. 5. Patterns formed by multiple particle systems when 
both attraction and repulsion are important. All images are 
acquired at f=20 Hz. (a) 3 particles at T — 3.0. (b) 7 par- 
ticles in a stable hexagon at V — 3.0. (c) 7 particles in sta- 
ble time-periodic motion at F = 3.7. (d) Many particles at 
r = 3.7, where the central particles are in contact and exterior 
particles are held apart by the hydrodynamic repulsion, (e) 
At r = 3.9 all the particles have separated and form a bound 
liquid, (f) T = 5.3, apparently chaotic. The time-dependent 
character of states (c) through (f) are revealed in on-line an- 
imations Jig ]. 

accelerations. Here the onset of a non-zero separation be- 
tween particles occurs at higher acceleration than for the 
two particle case. As shown in Fig. 5(d) (r = 3.7), par- 
ticles near the center of the cluster can remain in contact 
even when particles near the periphery have separated. 
One reason for this is that the vertical motion of particles 
is significantly smaller at the center of the cluster than at 
the edges, causing weaker repulsion. At T = 3.9, (Fig 5e) 
the inner particles have separated and the system forms 
a bound state with weakly chaotic motion of all the par- 
ticles, a "mesoscopic liquid". One interesting feature of 
this state is that the preferred distance between particles 
results in the formation of closed shells. At the param- 
eters of Fig 5(e) there are two more particles than fit 
in the closed shells, and these "valence" particles move 
more freely than the others. Further increase of T results 
in increased interparticle distance and more rapid parti- 
cle motion (Fig 5f). Animations of these time-dependent 



states are available [ p"5| |. 

We emphasize that the chaotic motion of particles 
in the clusters is not a result of external effects such 
as roughness of the driving surface. A single particle 
bounces regularly with possibly a slow drift due to im- 
perfect leveling or slightly non-vertical vibrator motion. 
Multiple particles exhibit complicated motion as a result 
of the nonlinearity of the flow of the interstitial fluid. 

In summary, we have shown that fluid mediated in- 
teractions between particles can be delicately tuned to 
produce a great variety of ordered and time-dependent 
dynamical states of these many body systems. 
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